Pattern formation in colloidal explosions* 
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We study the non-equilibrium pattern formation that emerges when magnetically repelling col- 
loids, trapped by optical tweezers, are abruptly released, forming colloidal explosions. For multiple 
colloids in a single trap we observe a pattern of expanding concentric rings. For colloids individually 
trapped in a line, we observe explosions with a zigzag pattern that persists even when magnetic 
interactions are much weaker than those that break the linear symmetry in equilibrium. Theory 
and computer simulations quantitatively describe these phenomena both in and out of equilibrium. 
An analysis of the mode spectrum allows us to accurately quantify the non-harmonic nature of the 
optical traps. Colloidal explosions provide a new way to generate well-characterized non-equilibrium 
behaviour in colloidal systems. 

PACS numbers: 82.70.Dd, 47.54.-r 



Introduction. - Pattern formation is an important 
and widespread phenomenon in the natural world and 
provides striking examples of order produced by non- 
equilibrium processes [1]. For instance, the confinement 
of charged particles to spatially localized traps can lead to 
one-dimensional Wigncr crystals or the spontaneous for- 
mation of zigzag or helical particle patterns. Such con- 
fined particle geometries are often used for ion storage 
and quantum-computing experiments [2, 3]. But when 
the confinement is too weak, the particles can escape 
from the traps. Often, particle trajectories produced 
from such "explosions" proceed in an incoherent manner, 
but in the case of long-ranged interactions, more inter- 
esting spatio-temporal patterns may form. In this letter, 
we investigate the dynamical behaviour of strongly in- 
teracting superparamagnetic colloidal particles that have 
first been confined with an optical laser trap. When the 
trap is abruptly turned off, the long-ranged repulsive 
inter-colloidal interactions generate a "colloidal explo- 
sion" with characteristic non-equilibrium patterns that 
depend on the initial conditions and on the strength of 
the interactions. 

Colloidal systems offer the unique ability to simultane- 
ously visualize and carefully control the non-equilibrium 
behaviour using external fields, which allows for detailed 
comparisons between experiments and theory [4] . Recent 
examples include non-linear instabilities in sedimenting 
suspensions [5], dynamic lane formation in oppositely 
charged particles under electric fields [6], driven dislo- 



cation nucleation [7] and stochastic resonance [8]. Here 
we exploit the ability to exquisitely tune inter-colloidal 
magnetic interactions through the application of an ex- 
ternal magnetic field [9] while simultaneously placing col- 
loids into a well-defined initial configuration using op- 
tical tweezers. The single-particle trajectories are di- 
rectly monitored by video- microscopy [10]. To explore 
the physics of colloidal explosions we study two basic 
geometries, namely the one-dimensional (Id) and two- 
dimensional (2d) configurations shown in Fig. 1. Even 
thought these geometries are fairly simple, they exhibit 
non-trivial non-equilibrium behaviour. This study opens 
up the possibility of exploring many other geometries and 
explosion patterns. 

Experimental and simulation details. - We use 

superparamagnetic latex spheres of radius a = 1.35 /Ltm 
(Dynabeads, Invitrogen) in a water solvent, contained in 
a 200 /urn thick quartz glass sample cell. The gravita- 
tional length of the particles is much smaller than their 
diameter so after sedimentation this is effectively a 2d 
system. The application of a perpendicular magnetic 
field induces long-range repulsive magnetic interactions 
between the colloids of the form 
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where \x§ is the permeability of free space, B is the mag- 
netic field, and r mn = |r m — r„| is the distance between 
particles n and m, and the colloids have a magnetic sus- 
ceptibility of x= 3.95xl0~ 12 A m 2 /T [11]. The particle 
trajectories are obtained using video- microscopy [10, 12]. 

The colloids are trapped using acousto-optical- 
deflection controlled optical laser tweezers (see e.g. [8]). 
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Figure 1: (Color online) (a) Microscopy images (92 x 23 /xm 2 ) showing an exploding Id array of N = 19 magnetic particles 
with a spacing L = 4 pm upon removing the optical traps at t — s at an external magnetic field B = 1.90 mT and stiffness 
ko = 0.37 ± 0.01 pN/pm. From top to bottom: t < s, 0.2s, Is, and 5s. (b) The experimental particle trajectories of the 
exploding chain (top) and for an exploding 2d colloidal disc (bottom) of JV = 28 particles at a magnetic field of B = 2.25 mT. 
The color code indicates the time in seconds, (c) Brownian dynamics trajectories for an exploding chain (top) and for a colloidal 
disc compare well to the experiments. 



The trapping potential for particle n is modeled as a 



V* rap = Vo 



1 — exp 
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(2) 



which takes into account, in a generic way, the fact that 
the trap has a finite range of attraction. Here, r„ is the 
particle position, R„ = nLx is the trap position, 5r n = 
r„ — R„ and Vq is the depth of the potential well. The 
softness of the trap is characterized by the dimensionlcss 
parameter 



koL 2 



(3) 



When a = the trapping is purely quadratic, Vjj rap = 
koSr^/2, for all Sr n , whereas for a > the trap is 
quadratic for small 5r n , but becomes increasingly non- 
harmonic at larger Sr n and has a finite height Vq above 
which the particle can escape from the trap. It is impor- 
tant to note that ko and Vq can be changed in the ex- 
periment but that their ratio remains virtually constant 
and therefore characterizes the optical trap [13, 14]. 

The maximum particle velocities observed in the ex- 
periments were of order 2 /im/s in the explosion so that 
the Reynolds number Re ~ 10~ 5 is small. Simple over- 
damped Brownian Dynamics (BD) simulations with a 
single fit parameter to the field B in Eq. (1) were found 
to closely mimic the experimental trajectories as demon- 
strated in Fig. 1(b) and (c). Combining the repulsive and 
trapping potentials, Eqs. (1) and (2), at finite tempera- 
ture T leads to the following dimensionless equation of 
motion 

in = 3b 2 ^ - 5r « CX P (-f 6r l) +^n- (4) 



Here and afterwards, the length, the energy, and the time 
are expressed in terms of of L, kgL 2 , and 6irr/a/ko (r/ is 
the dynamic viscosity of the solvent), respectively. The 
stochastic force s(t) obeys the properties (s n (i)) = 
and (s n (t)s n i(t')) = 5 nn >5(t — t') with the first S being 
Kronecker's delta function and the second being Dirac's. 
Apart from the parameter characterizing the softness of 
the trapping potential, a, the dynamics is governed by 
dimensionless parameters 



inknL 5 ' 



k B T 
knL 2 ' 



(5) 



These two additional parameters describe respectively 
the intensity of the magnetic field and the intensity of 
thermal motion. The fact that we work at a given tem- 
perature T, and at certain values of L and ko, fixes the 
value of e. To integrate Eq. (4), we apply a standard 
algorithm [15]. 

Explosion of a disc. - We first describe the 2d ge- 
ometry dealing with explosions from a disc composed of a 
variable number of colloids, as shown for an example with 
N = 28 particles in Fig. 1(b). To create such a disc con- 
figuration in our simulations, all the colloids are captured 
by a wide single trap, R„ = (0, 0). However, as we argue 
below, our experimental optical trap is characterized by 
the value a « 30, implying that the effective entrapment 
range is such that the trap cannot easily hold more than 
one colloid at a time. In the experiments, we therefore 
gathered all the colloids together into a disc, which was 
easily done with a single trap in the absence of the mag- 
netic field. As soon as the colloids have been gathered 
in the disc, the trap was switched off and, at the same 
time, the magnetic field was quickly switched on. Typ- 
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Figure 2: (Color online) Particle trajectories for explosions 
from two-dimensional discs of N = 28 particles with and with- 
out delay. The values of parameters b — 0.274 and a — are 
chosen to provide the best fit to the experiment, (a): no de- 
lay, At = 0; (b): delay time At ~ td/10; (c): delay time 
At « to/5. 



ically, this procedure is done within the time At < 1 s. 
As the characteristic diffusion time for the particles used 
in experiments is td ~ a 2 / D rj 10 s (where D is the 
particle diffusivity) , the condition At < Trj ensures that 
the gathered colloids have not significantly diffused away 
from their initial configuration. 

We observe a clear pattern of concentric rings as the 
colloids move to minimize the magnetic repulsion be- 
tween them (see Fig. 1(b) and (c)). The velocities v 
of particles are maximal in the beginning of the explo- 
sion and decay as the particles move apart. As the re- 
pulsive interactions are long ranged, in an unbounded 
domain without thermal fluctuations this process never 
stops and the particles move out to infinity. Neverthe- 
less, at finite temperature a characteristic explosion time 
T e can be defined as the time to reach a regime of motion 
with a small enough velocity such that the Peclet number 
Pc = av/D ~ 1. At smaller velocities, diffusion starts to 
dominate. 

In Fig. 2 we illustrate simulated trajectories for explo- 
sions of a disc of N = 28 particles for different delays 
Ar, in which the explosion time r e ~ tq/6. Note that 
the trapped colloids can come very close to each other. 
To ensure that the colloids do not overlap in simulations, 
we have additionally included in Eq. (4) steep repulsive 
interactions of the Weeks-Chandler- Andersen form [16]. 
These results confirm that our way of creating a disc in 
the experiment is practically similar to the use of a sin- 
gle wider trap in the simulations as described above. We 
emphasize that the results remain robust even for delay 
times that are longer than in the experiment and com- 
parable with the time of explosion, as, e.g., in Fig. 2(c). 

In figure 3 we show simulation trajectories for discs 
composed of different number of colloids. The initial 
configuration in the form of a 2d disc can be considered 
as a pattern of concentric rings that helps explain the 
non-equilibrium pattern seen in the explosions. As N in- 
creases, the number of rings grows, with a single particle 
remaining near the centre for N = 6, marking the begin- 
ning of a second ring, and then again at N = 16, marking 
the beginning of a third ring, and again at N = 31, mark- 



ing the beginning of a fourth ring. The pattern formation 
is caused by the initial shell-like ordering of the particles 
that is driven by the combination of repulsive magnetic 
fields and confinement [9]. The particles in each shell 
move outwards in the same manner until the repulsive 
magnetic interactions between the particles are less than 
ksT {ks is Boltzmann's konstant) and diffusion begins 
to dominate. We note that this colloidal system is rem- 
iniscent of Coulomb explosions induced by strong laser 
fields stripping off the electrons of molecules and atomic 
clusters [17, 18]. 

Explosion of a chain. - Whereas the broad fea- 
tures of the colloidal disc explosions are fairly straight- 
forward to explain, the behaviour of the Id configura- 
tion is more subtle. This is not unexpected since Id 
chains of interacting particles form some of the rich- 
est and most important models in theoretical physics. 
Examples include the Frenkcl-Kontorova model [19], Id 
models of Wigncr [20, 21], colloidal [22, 23] and microflu- 
idic droplets [24, 25] crystals and other closely related 
systems such as polymer chains [26] and dusty plasmas 
[27]. 

The experiment shown in Fig. 1(a) and (b) is for N = 
19 colloids placed with a lattice spacing of L = 4.00 /im 
apart. When the optical traps are removed in the pres- 
ence of the magnetic field, the particle chain explodes and 
the evolving pattern shows an almost perfect zigzag sym- 
metry. To analyze the underlying equilibrium behaviour 
in the presence of a magnetic field, but before the optical 
traps are switched off, i.e. before the explosion, we begin 
by treating the deterministic (e = 0) case for N — > oo 
colloids. For a strong enough inter-colloidal repulsions, 
and for strong enough optical traps, we expect that the 
linear symmetry will be broken, leading to an equilibrium 
transition from a linear to a quasi-onc-dimensional zigzag 
state as, e.g., in the Frenkel-Kontorova model [19], dusty 
plasmas [27], quantum wires [28] and Id Wigncr crys- 
tals [21]. The zig-zag state can be described by the order 
parameter h, where h n /L = nx + ( — l) n (/i/2)y is the 
equilibrium displacement from the centre of the trap of 
particle n. Combining the induced interaction (1) with 
the laser traps potential (2) leads to an interaction energy 
(per particle) U(h) that can be written in the dimension- 
less form 

U(h)=^±f m (h) + 1 - e ^- ah2/8 \ (6) 
Z — ' a 

m— 1 

with f m (h) = (m 2 + Pm h 2 )- 1 / 2 and p m = [1 - (-l) m ]/2. 
The equilibrium displacement /i* can be found as a func- 
tion of b and a by minimizing the potential (6). The 
critical field b c at which the transition from the line to 
the zigzag state with h = h* ^ occurs is found to be 
b c = v/8/[93C(5)] w 0.288, with ((x) = £~ =1 mT x the 
Ricmann Zeta function. Close to the critical point, where 
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Figure 3: (Color online) Particle trajectories for explosion from two-dimensional discs composed of different number N of 
colloids. Parameters are b = 0.274, a = 0. 



\b — b c \ <C b c and ft, Cl, we obtain a square-root law 

^ ± jp| (7, 
y 6^(cv c - a) 

with a c = 635C(7)/[31C(5)] w 19.9. We note that the 
nearest neighbor (NN) approximation, where only one 
term with m = 1 in relation (6) is retained, works very 
well: We obtain bf w = l/\/l2 « 0.289 and af w = 20. 

We note that the existence of the equilibrium zig-zag 
state is determined by the value of the softness param- 
eter, a. For potentials with a < a c , the Id line state 
is stable for b < b c , whereas for b > b c the zig-zag 
state is stable and the displacement h* grows with in- 
creasing b. However, for < a < a c , increasing the 
field further eventually leads to an instability at a larger 
field 6»* where the equilibrium zig-zag state is unstable, 
i.e. the barrier to the particles leaving the wells disap- 
pears. Within the NN approximation we find that = 
bc N (u? N A*) * exp[- (af N - a)/16], although finite tem- 
perature means particles can escape at lower fields, 
K* = Ml - (l/40)exp[(a c - a)/12](30k B T/V ) 2 ^}. 
For softer potentials with a > a c , there is no transition 
to an equilibrium zigzag state. In this case, the parti- 
cles experience a colloidal explosion and leave their wells 
before an equilibrium zigzag state can occur. In the ab- 
sence of thermal fluctuations, increasing the field would 
yield an explosion of the Id line state at b = b c . At fi- 
nite temperature, particles can escape at smaller fields, 
K « b c [l - y/{l-a c /a)k B T/2Vo]. 

The question then arises whether for the experiments 
set-up used in Fig. 1, a is small enough to allow an 
equilibrium zig-zag transition. To test this we measured 
the probability density functions for displacements along 
(P(Sx)) and perpendicular to (P(5y)) the trapped chain 
for different magnetic fields. As can be seen in Fig. 4, 
P{Sy) widens for increasing magnetic fields, but always 
retains its single-peaked structure. This confirms that 
the system is not in an equilibrium zigzag state, despite 
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Figure 4: (Color online) Probability density distribution in 
the transversal direction, P(5y), for increasing magnetic field 
B. The inset shows respectively the standard deviations a x 
(solid line) and a y (dashed line) of P(8x) and P(Sy) as func- 
tions of magnetic field. Points represent experimentally mea- 
sured data, the lines provide the fit obtained via BD simula- 
tions with ft = 30. 

the zigzag symmetry in the explosion pattern originating 
from these states. The inset of Fig. 4 presents respec- 
tively the standard deviations a x and a v for distributions 
P(Sx) and P(Sy) as functions of the magnetic field. The 
observation that a x decreases and a y increases upon in- 
creasing the magnetic field points to respectively harden- 
ing of the longitudinal and softening of transverse normal 
modes [21]. 

We also measured the phonon dispersion relations by 
tracking the particle displacements u(h„) from their 
equilibrium positions h„. The Fourier transforms u(q) 
of the displacement vectors are directly related to the 
dynamical matrix D^ u (q) [10]: 

(u* fl (q)u v (q)) = k B TD-l(q), (8) 

where the average is over all independent configurations. 
The eigenvalues of D^ v (q) yield the normal mode spring 
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Figure 5: (Color online) The longitudinal and transverse nor- 
mal mode spring constants fc"' x for different magnetic fields 
as a function of the wave number q. Points are experimental 
data. Lines are from BD simulations with a = 30. 



constants. Figure 5 shows the phonon-dispersion rela- 
tions as a function of q for different values of the magnetic 
field and the hardening and softening of the longitudinal 
and transverse modes respectively is indeed evident. 

We analyzed these results further using BD simula- 
tions. Given the spring constant, the interparticlc dis- 
tance L and the temperature, there are only two param- 
eters left to determine, the ratio between B and b and the 
potential softness parameter a. A fit to the data in both 
Fig. 4 and Fig. 5 yields a single proportionality constant 
B/b and also fixes a = 30 ± 5. There is good quantita- 
tive agreement with the data in both figures. The values 
of b < b c , suggest that we are below the critical field, 
and furthermore the fact that a > a c suggests that for 
this experimental system no equilibrium zigzag state is 
possible. 

We note that varying a affects the transverse and lon- 
gitudinal modes differently. For a = 0, A;"' x (g = 0) = 1 
(6 < b c ), whereas for a ^ we obtain at q = 0: 
fc" =/= k 1 - ^ 1, in agreement with experiment. Note that 
k^' 1 - are in units of fcrj. These deviations can be traced 
back to the softness of the trapping potentials. At finite 
temperature, the particles explore the non-harmonicity 
of the potential, leading to what appears to be a weaker 
effective spring force. As before, the accurate NN ap- 
proximation helps illuminate these results. The normal 
modes correspond to u(h„) oc nexp(— At + inq), where 
A is the decay rate, q G [0, tt] is the wave number, and 
n = x or y. For the line state we obtain, in dimensionless 
form, 
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from which we draw three important conclusions. 



Figure 6: (Color online) Particle trajectories showing zig-zag 
explosions for chains with different number of colloids: iV = 9 
(a), N = 19 (b), and N = 29 (c). d) A pattern with defects 
for iV = 29. Parameters are: b = 0.86 c , a — 30. 



First, this result shows hardening and softening of the 
spring constants with the field in respectively the longi- 
tudinal and transverse directions. Since k"' X oc &~y, this 
is consistent with the inset of Fig. 4. Note that the de- 
pendencies of fell- 1 on q at a given b provided by Eqs. (9) 
and (10) are in qualitative agreement with those obtained 
via relation (8) and presented in Fig. 5. Second, Eqs. (9) 
and (10) show that as b goes beyond 6 C , the mode that 
first becomes unstable is the zigzag mode, q = tt. Third, 
in a system with a finite number of particles the spec- 
trum of decay rates is discrete. There are other eigenval- 
ues present as well (cf. markers in Fig. 5), so that the 
zigzag mode is better separated from other eigenmodes 
in chains with smaller number of beads, N. Hence, while 
the perfect zigzag pattern is the most probable in explo- 
sion of shorter chains, the probability of imperfections in 
this pattern increases with the chain length. We indeed 
observed this effect in simulations of longer chains. 

In Fig. 6 we show explosion patterns for Id chains of 
different lengths. We fix the magnetic field to a subcriti- 
cal value b = 0.8b c and the softness parameter to a = 30 
(as in the experiments) which guarantees that the ex- 
plosion occurs from the equilibrium line state. As can 
be seen from the trajectories, similar explosion patterns 
with the zig-zag symmetry can be found in shorter and 
longer chains. However, for the longer chains, explosion 
patterns with defects also occur, see Fig. 6(d). The sim- 
plest defect corresponds to a pair of neighbouring colloids 
shooting out in the same direction. Defects appear be- 
cause phonon modes other than the zigzag mode play an 
increasingly important role with increasing N. 

For systems with a < a c an equilibrium zigzag state 
can always be induced when > b > b c . In contrast to 
systems without an equilibrium zigzag state, we observe 
perfect zigzag symmetry in the explosions for all we 
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investigated. The different behaviour results because in 
this case the non-equilibrium pattern simply reflects the 
broken symmetry in equilibrium. 

The fact that a single fit to B/b and a provides an 
accurate fit to all the data from Figs. 4 and 5, combined 
with the single-peaked probability distributions, strongly 
suggests that the explosions observed in Fig. 1 are for 
the situation where there is no equilibrium zigzag pat- 
tern. The zigzag pattern we observe is entirely due to 
a non-equilibrium pattern formation effect that shadows 
the hidden equilibrium symmetry breaking. 

Finally, we also considered the role of hydrodynamic 
interactions (HI) by using 3d Stokesian dynamics. We 
find very similar behaviour to that seen with BD simula- 
tions that neglect the HI. The real HI will be modulated 
by the surface, but given the small overall effect of HI, 
we argue that explicitly including them is not important 
for the colloidal explosions we studied. 

Conclusions. - In summary, we have exploited the 
ability to carefully control and characterize a colloidal 
model system to generate non-equilibrium pattern for- 
mation in colloidal explosions. For a 2d initial geome- 
try, we observe patterns with expanding concentric rings 
that can be explained by the initial shell-like ordering. 
For the Id geometry, we observe explosions with a zigzag 
pattern at fields strengths well below those that would 
break the linear symmetry in equilibrium. Furthermore, 
a quantitative comparison to phonon dispersion relation- 
ships allows us to characterize the non-harmonic nature 
of the traps. 

More generally, we have introduced "colloidal explo- 
sions," a new way to generate well-characterised non- 
equilibrium behaviour in colloidal systems. It should be 
possible to create such explosions with a range of different 
kinds of confining potentials and repulsive inter-particle 
interactions. This methodology can be applied to a wide 
variety of other geometries, leading to potential applica- 
tions, for example, in microfluidics. 
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